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ABSTRACT 

Assuming an adiabatic evolution of a Gamma-Ray Burst (GRB) remnant interacting with an external medium, 
we calculate the injection, cooling, and absorption break frequencies, and the afterglow flux for plausible orderings 
of the break and observing frequencies. The analytical calculations are restricted to a relativistic remnant and, in 
the case of collimated ejecta, to the phase where there is an insignificant lateral expansion. Results are given for 
both a homogeneous external medium and for a wind ejected by the GRB progenitor. 

We compare the afterglow emission at different observing frequencies, for each type of external medium. It is 
found that observations at sub-millimeter frequencies during the first day provide the best way of discriminating 
between the two models. By taking into account the effect of inverse Compton (IC) scatterings on the electron 
cooling, a new possible time-dependence of the cooling break is identified. The signature of the up-scattering 
losses could be seen in the optical synchrotron emission from a GRB remnant interacting with a pre-ejected wind, 
as a temporary mild flattening of the afterglow decay. The up-scattered radiation itself should be detected in the 
soft X-ray emission from GRB remnants running into denser external media, starting few hours after the main 
event. 

Subject headings: gamma-rays: bursts - methods: analytical - radiation mechanisms: non-thermal 



1. INTRODUCTION 

One of the most important issues regarding Gamma-Ray 
Bursts (GRBs) is the nature of the object that releases the rela- 
tivistic ejecta generating the high energy emission of the main 
event and the lower frequency emission during the ensuing af- 
terglow. Some insight about the GRB progenitor can be ob- 
tained from the properties of the circum-burst medium, which 
can be inferred from the features of the afterglow emission. If 
the ejecta is expelled during the merging of two compact ob- 
jects (Meszaros & Rees 1997b), it is expected that the medium 
surrounding the GRB source is homogeneous. However, if a 
collapsing massive star (Woosley 1993, Paczyhski 1998) is the 
origin of the relativistic fireball, the circum-burst medium is the 
wind ejected by the star prior to its collapse, whose density de- 
creases outwards. The two models differ in the dependence on 
radius of the particle density of the circum-burst medium which 
the GRB remnant interacts with, and in the value of this den- 
sity at the deceleration length-scale. The former modifies the 
rate of decline of the afterglow, while the latter determines the 
overall afterglow brightness. Therefore it is possible to corre- 
late afterglow emission features to a specific type of external 
medium. 

Significant work in this direction has been done by many 
researchers. The two afterglows that exhibited breaks consis- 
tent with the effects arising from strong collimation of ejecta - 
GRB 990123 (Kulkarni et al. 1999a) and GRB 990510 (Stanek 
et al. 1999, Harrison et al. 1999) - indicate that the external 
gas was homogeneous (recent work by Kumar & Panaitescu 
2000 shows that jets interacting with winds cannot produce 
sharp breaks in the afterglow light-curve). The optical emis- 
sion of three afterglows had a steeper than usual decline. GRB 
970228 decayed as T~ 17 after the subtraction of an underlying 
supernova emission (Reichart 1999, Galama et al. 2000). The 



light-curve of GRB 980326 fell off as T~ 2A (Groot et al. 1998) 
and an emission in excess of the early time extrapolation was 
detected ~ 20 days after the main event, indicating a super- 
nova contribution (Bloom et al. 1999). A T~ 2 decay was ob- 
served for the afterglow of GRB 980519 (Halpern et al. 1999). 
Such steep declines can be produced either by a fireball inter- 
acting with a pre-ejected wind (Chevalier & Li 1999) and an 
electron index around 3, or by a narrow jet expanding later- 
ally in a homogeneous external medium and an electron index 
slightly larger than 2. Chevalier & Li (1999) found that the ra- 
dio emission of the afterglow of GRB 9805 19 is consistent with 
an external wind; however Frail et al. (2000) point out that the 
interstellar scintillation present in the radio data does not allow 
ruling out the jet model. Nevertheless, the existence of super- 
novae associated with GRB 970228 and GRB 980326 points 
toward a massive star as the origin of these bursts, implying 
a pre-ejected wind as the external medium. From the analysis 
of the optical radio emission of the afterglow of GRB 970508, 
Chevalier & Li (2000) conclude that the surrounding medium 
was a wind. Frail, Waxman & Kulkarni (2000) argue that the 
same radio afterglow can be explained by a homogeneous ex- 
ternal medium. 

In this work we investigate the differences between the light- 
curves of afterglows arising for each type of external medium, 
with the aim of finding ways for distinguishing between the 
two models. This study is done within the usual framework of 
a relativistic remnant interacting with a cold external gas. As 
the fireball is decelerated, a shock front sweeps up the external 
gas, accelerating relativistic electrons and generating a mag- 
netic field in the shocked gas. We ignore the emission from 
electrons accelerated by the reverse shock which propagates 
through the ejecta at very early times. At optical wavelengths 
this emission is short lived, lasting up to few tens of seconds af- 
ter the main event (Sari & Piran 1999), but it could be important 
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for the radio emission until few days (Kulkarni et al. 1999b). 

Analytical afterglow light-curves for spherical remnants in- 
teracting with homogeneous external media have been previ- 
ously published by Sari, Piran & Narayan (1998). Features of 
afterglows from spherical fireballs, such as peak flux, break fre- 
quencies, and time evolution of fluxes at a fixed frequency, have 
been studied by Meszaros & Rees (1997a), Waxman (1997b), 
Wijers & Galama (1999), and Dai & Lu (2000) for homoge- 
neous media, and by Chevalier & Li (2000) for pre-ejected 
winds. In this work we present and compare analytical and 
numerical light-curves at various observing frequencies, cover- 
ing all cases of interest, for both types of external media, tak- 
ing into account the differential arrival-time delay and Doppler 
boosting due to the spherical shape of the source. We take into 
account first order IC scattering, calculate its effect on the elec- 
tron cooling and on the afterglow synchrotron emission, and 
study briefly the high-energy emission resulting from the up- 
scattering of synchrotron photons. The possible importance of 
IC scatterings for the early afterglow emission was pointed out 
by Waxman (1997a) and Wei & Lu (1998). 



2. SIMPLE DYNAMICS OF RELATIVISTIC REMNANTS 

For the calculation of the afterglow emission it is necessary to 
know how the remnant Lorentz factor V evolves with observer 
time T, as all other quantities that appear in the expression of 
the spectral flux are functions of T and of the remnant radius 
r and external medium density n(r). We shall assume that the 
remnant is adiabatic, i.e. the energy carried away by the emit- 
ted photons is a negligible fraction of the total energy of the 
remnant. This assumption is correct if the energy density of the 
electrons accelerated at the shock front is a fraction e e <C 1 
of the total energy density in post-shock fluid or if most of the 
electrons are adiabatic, i.e. their radiative cooling timescale ex- 
ceeds that of the adiabatic losses due to the remnant expansion. 

Assuming that the internal energy of the ejecta is negligible 
compared to its rest-mass energy and that the ratio internal-to- 
rest mass energy in the energized external medium is T — 1 
(i.e. its "temperature" tracks that of the freshly shocked gas), 
conservation of energy leads to 



m(r)T 2 + M fb T - [m(r) + M fb T ] = 



(1) 



where Mfb and Tq are the initial mass and Lorentz factor of the 
fireball (whose energy is E — MfbTo) and 



m(r) 
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m p n(r)r 3 



(2) 



is the mass of swept-up material (m p being the proton's mass). 
The external medium particle density is 



n(r) — Ar s , 



(3) 



with s = for a homogeneous medium and s = 2 for a wind 
ejected by the GRB progenitor at a constant speed. 

Equation (Q) is valid if the remnant is spherical, but can 
also be used for collimated ejecta when the lateral spreading 
(Rhoads 1999) is insignificant if the quantity E above is de- 
fined as the energy the fireball would have if it were spherical. 
Throughout this work we shall assume that the remnant is a jet 
with an initial half-angle larger than > 20°, in which case the 
sideways expansion is negligible during the relativistic phase. 



The following analytical calculations of the afterglow emission 
can be extended to sideways expanding jets and non-relativistic 
remnants by first determining T(r). The set of coupled differen- 
tial equations describing the evolution of the jet Lorentz factor 
and its opening can be solved analytically for s = (Rhoads 
1999). The lack of a good approximation for the jet dynam- 
ics in the case of pre-ejected winds is the main motivation for 
restricting the following analytical calculations to spherical or 
wide-angle remnants. 

The solution of equation (pi) is 



r(r) 



1 + (2a: 3 -vr ) 2 - 1 
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where x is the radial coordinate r scaled to 



ro 



E 



Air m v c 2 AT?. 



l/(3-s) 



(4) 



(5) 



the deceleration length-scale, at which m(ro) = E/(c 2 Tq) = 
Mfh/T . The result given in equation (Q) is also valid in 
the non-relativistic regime. For x <C 1 T < To, while 
for 1 < x < x nr we find T = a;~ (3 ~ s)/2 r . Here 
Xnr = (r 2 l /3) 1 / ( - 3 ~ s " ) marks the end of the relativistic regime: 
r(x nr ) = 2. For the ease of analytical calculations we shall 
assume that the power-law behavior of T lasts from x = 1 to 

The Lorentz factor given in equation (Q) represents a "dy- 
namical" average of the Lorentz factors at which different re- 
gions of the shocked remnant move (the Blandford - McKee 
solution). The Lorentz factor of the shock front that propagates 
into the external gas r s ;, = \[2 T, with T given by equation 
(Q), matches that given in equation (69) of Blandford & McKee 
(1976) for the power-law regime 1 <C x <C x nr if E is multi- 
plied by (17 - 4s)/(12 - 4s). This correction factor (17/8 for 
s = and 9/4 for s — 2) will be used in the following results. 

The constant A in equation (|j) is the number density n» of 
the external homogeneous medium for s = 0, while for s = 2 

A = = 3.0 x 10 35 A* cm- 1 , (6) 

where M is the mass loss rate of the massive star that ejected 
the wind at constant speed v, and A was scaled to 



A± 



M/lO^Meyr- 1 
v/10 3 kms- 1 



(7) 



as in Chevalier & Li (2000) for a Wolf-Rayet star. For fur- 
ther calculations it is convenient to use equation (||) in the form 

n(r) = no(r /ro)~ s where no = rc* for s = 0, while for s = 2 

n = 1.9 x 10 4 E^Tl 2 Al cm" 3 (8) 
is the wind particle density at the deceleration radius (eq. [Q|): 

(s = 2) r Q = 4.0 x 10 15 E 53 T 2 2 A^ cm . (9) 

The usual notation C„ = 10~™C is used throughout this work. 

Note that, for the reference values used here, the deceleration 
radius in the wind model is 1 .5 orders of magnitude lower than 
that for a homogeneous external medium: 

(s = 0) r Q = 1.3 x 10 17 £4 /3 r o i /3 <o /3 cm . (10) 
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and even smaller for higher initial fireball Lorentz factors or 
slower winds (i.e. a larger parameter A*). If GRBs are due to 
internal shocks occurring in unstable relativistic fireballs (Rees 
& Meszaros 1994, Paczynski & Xu 1994, Piran 1999), then the 
external shock resulting from the interaction of the fireball with 
the pre-ejected (non-relativistic) wind may occur before the in- 
ternal shocks are over. In this case successive internal collisions 
occur when faster parts of the ejecta catch up with the decelerat- 
ing leading edge of the fireball, a scenario suggested within the 
s = model by Fenimore & Ramirez-Ruiz (2000), but which is 
more likely to happen if the external medium is the gas ejected 
by a massive star. The GRB itself would then exhibit the erratic 
variability characteristic of internal shocks until a time of the 
order 

^-we^t^a; 1 *> (ii) 

Cl o 

(which has a strong dependence on T ), after which there may 
be significant emission from internal shocks on the outermost 
part of the fireball and from the external shock that plows 
through the external gas. The former mechanism generates 
pulses of increasing duration as the fireball expands, while the 
later leads to a continuous emission. 

The time T when the observer receives a photon emitted 
along the line of sight toward the fireball center can be calcu- 
lated by integrating 



dT = (1 - P)dt 



1 dt 
2P 



(12) 



where (3 is the shocked fluid speed and t = r/c is the time 
measured in the laboratory frame. Approximating the solu- 
tion given in equation (0) with T = Tq for x < 1 and 

T = x~( 3 ~ s ^ 2 To for 1 < x < x nr , one obtains 



T = T (x 4 



To 



'"o 



(13) 



2(4 -a) cTg • 

Note that T given in equation ( |l3| ) is the earliest time a photon 
emitted by the remnant at time t can reach the observer. Pho- 
tons emitted by the fluid moving at an angle 9 = 1/T off the 
center-observe axis arrive at T = t(l — cos 9) = t/2T 2 , which 
is a factor 4 — s larger than the time corresponding to 9 = 0: 
T = i/2(4- s)T 2 . 

From equation ( pj| ) r(T) can be found and then substituted 
in the expressions for T(r) and n(r) to obtain these quantities 
as a function of the observer time. For the power-law phase the 
results are 



( S = 0) 



T{T)=&.ZEli & n-l' & Tf /& 



3. BREAK FREQUENCIES 

Within the synchrotron emission model there are three ex- 
pected breaks in the afterglow spectrum: (i) an injection break, 
at the synchrotron frequency vi at which the bulk of the elec- 
trons injected by the shock front radiate, (ii) a cooling break, 
at the synchrotron frequency v c of electrons whose radiative 
cooling time equals the expansion timescale, and (Hi) an ab- 
sorption break, at v a below which the synchrotron photons are 
absorbed by electrons in free-free transitions in a magnetic field 
(synchrotron self-absorption). 

The break frequencies can be calculated if the distribution of 
the injected electrons and the strength of the magnetic field are 
known. The distribution of the injected electrons is assumed 
to be ATi(j) oc 7~ p starting from a minimum random Lorentz 
factor given by 



7i = I% e (r-i), 



(19) 



where m e is the electron mass. The energy carried by this elec- 
tron distribution is a fraction ^z^£ e of the total internal energy. 
The post-shock magnetic field strength in the co-moving frame 
is given by 



B 



— = £Bm p c 2 n' e (T-l) = Ae B m p c 2 n(r)(T-l) \T + - 

(20) 

where eb is the fractional energy carried by the magnetic field 
and n' e is the co-moving frame electron density behind the 
shock front. Equations ([19]) and ( [20| ) are based on that the in- 
ternal to rest-mass energy density ratio in the shocked fluid is 
r — 1; the derivation of the latter equation also used that the co- 
moving particle density is AT + 3 times larger than that ahead 
of the shock. 



3.1. Injection Break 

Using the relativistic Doppler factor 2T corresponding to the 
motion of the source toward the observer (i.e. 9 = 0), the syn- 
chrotron emission from a power-law distribution of electrons 
peaks at the observer frame frequency 



(14) 



2ir m e c 



-y 2 BT = 8.4 x 10 6 x p -f 2 BT Hz 



(21) 



(s = 0) r(T) = 8.2 x 10 17 E^n'^T^ cm , (15) 
n = (constant) for a homogeneous external medium and 



(s = 2) T(T) = 7.9 2Sg£ A* TJ" , (16) 

(s = 2) r(T) = 6.4 x 10 17 El[ 2 A^ 1/2 T^ 2 cm , (17) 



(s = 2) 



n(T) = 0.73 E^AlTZ 1 cm" 



(18) 



for an external wind, Td being the observer time measured in 
days. Note that, at least for the scaling values chosen here, T, 
r, and n have about the same values at T — 1 day in both mod- 
els. Also note that the above quantities (and thus the afterglow 
emission) are independent of the fireball initial Lorentz factor 

r . 



where the factor x v is calculated in Wijers & Galama (1999) for 
various values of the electron index p. We shall use x p — 0.52, 
which is strictly correct only for p — 2.5. With the aid of equa- 
tions @, @, @, @, and @ one obtains: 



(s = 0) Vi = 0.92 x 10 13 i4 /2 <-i4 /2 -2^ 3/2 Hz 



(22) 



(s = 2) Vi = 1.9 x 10 13 E^s%_ x E^_ 2 T d 612 Hz . 

(23) 

Note that Vi has the same scalings with the model parameters 
for s = and s = 2. The ratio of the two frequencies is \ r^. 
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3.2. Cooling Break 

The relativistic electrons cool radiatively through syn- 
chrotron emission and IC scatterings of the synchrotron pho- 
tons on a co-moving frame timescale 



t radii) = 



6n m P c 1 



Y + l Y + 1 CT e 7.B 2 ' 



(24) 



where t' sy (^) is the co-moving frame synchrotron cooling 
timescale of electrons of Lorentz factor 7, Y is the Compton 
parameter, and cr e the cross-section for electron scatterings[J 
Using equation (p0|), the Lorentz factor of the electrons that 
cool radiatively on a timescale equal to the remnant age 



cj r 5 



s cT 



can be written as 

3tt(5 - s) m e c 2 T 77(5 - s) 1 



7c 



y + l cr e 5 2 r (Y + l)e B nl>i 8 



(25) 



, (26) 



with rt in cm 3 . 

The observer-frame frequency v c of the cooling break is 



(s = 0) v c = 3.7xl0 14 E~ 11 ^(Y+iy'e-^Tr" Hz 

(27) 

(s = 2) i/ c = 3.5 x 10 14 Ell 2 A- 2 (Y+l)- 2 e-^lT¥ 2 Hz . 

(28) 

Hereafter we shall use the terminology "radiative electrons" 
for the case where the ji -electrons cool mostly through emis- 
sion of radiation (i.e. t' rad {~fi) < t' and 7 C < 7;), and we shall 
refer to "adiabatic electrons" if the 7;-electrons cool mostly adi- 
abatically (i.e. t' rad (ji) > t' and j t < j c ). 

3.2. 1 . Compton Parameter and Electron Distribution 

For the calculation of the Compton parameter Y, we take 
into account only one up-scattering of the synchrotron photons. 
Multiple IC scatterings of the same photon have an important 
effect on the electron cooling only if the Y param eter for sin gle 
scatterings is above unity. As shown in §3.2.2 and §3.2.3 this 
occurs if 1) eg £ 10~ 2 e e ,_i and if 2) T < 1 y , where T y is the 
time when Y falls below unity. 

Using equations ( |l9| ) and (^), it can be shown that 
IC scatterings of order higher than two are suppressed by 
the Klein-Nishina effect. For adiabatic electrons, a second 
IC scattering occurs in the Thomson regime if 3a) T <; 

1 E li Z <l' 9 z 2 °-\e 2 B 9 -i d for s = °< and if 3fo ) T ^ 
1 E^A* 1 ^ 4 £^-i£ 1 b 4 -4 d f°i" s = 2. If a second IC scattering 
is ignored in these cases then a higher energy component peak- 
ing around 1 GeV is left out, otherwise the synchrotron and first 
IC emissions remain unaltered if the electrons are adiabatic. 

The effect of second order up-scatterings is more important 
when ele ctron s are radiative, i.e. for 4) T < T r with T r calcu- 
lated in § 3.2.2, as in this case it reduces the intensity of the syn- 
chrotron and first stage IC components. With the aid of equa- 
tions (H) and © it can be shown that for s = a second 
up-scattering occurs in the Thomson regime and at T > 10 -2 d 



-2-3/2,7,-1/2 



if 5a) n > 25 E, 



7/11^-10/11 



£ B 2 !_^ cm 3 , while for s = 2 

1 The upjSca ttering of the Vi synchrotron photons on the 7 C - and 7^ -electrons 
in equation (E4I) is not reduced by the Klein-Nishina effect 



the second IC emission is not suppressed by the Klein-Nishina 
effect if 56) T > 0.05 E~ 3 1/3 Al^e^e 2 ^ d. 

Concluding, second stage up-scatterings can be ignored if the 
set of conditions 1), 2), 3) or 1), 4), 5) are not simultaneously 
satisfied. For a single up-scattering, the Compton parameter is 

Y=\l l 2 &r e = -T e [ 7 2 A/- e ( 7 )d7, (29) 

where Af e (7) is the normalized electron distribution and r e is 
the optical thickness to electron scattering, given by 



1 G e m(r) 1 



47r m p r 2 3 — s 



(30) 



If the injected 7^ -electrons cool faster than the timescale of 
their injection, then 7 C given by equation (^) is the typical elec- 
tron Lorentz factor in the remnant, and the electron distribution 
in the shocked fluid can be approximated by 



AC (r) oc s _ {p+1) 



7c < 7 < 7« 
It < 7 



(31) 



with p > 2. In the opposite case most electrons have a random 
Lorentz factor 7,, and the electron distribution is 



7 

K 1 -V-CP+ 1 ) 



7i < 7 < 7 C 



(32) 



3.2.2. Radiative Electrons 
For 7 C <g; 7; equations (El) and (pl|) lead to 



Y- = TflilcTe ■ 

Substituting 7 C with the aid of equation (Eq), one obtains 



(33) 



Y (Y 1) - 5 ~ 5 < m -° 2 ^ - 5 ~ s £e G4) 
r[r+ ' 8(3 -a) B 2 /8tt 8(3 - a) e B ' ( V 

where we used equations (|l9|) and (p0|). Therefore the Compton 
parameter during the electron radiative phase is 




2(3 -s) e B 



1 - 1 



(35) 



Hence the electron cooling is dominated by IC scatterings 
(i.e. Y r > 1) for eb < £~b, where 



16(3 - s) 



(36) 



Note that Y r is time-independent. Therefore the v c given by 
equations ( ]27| ) and ( ^8| ) decreases with time for s = and in- 
creases in the s = 2 model. Thus, for observations made at a 
fixed frequency, the electrons emitting at that frequency change 
their cooling regime from adiabatic to radiative in the case of a 
homogeneous external gas, and from radiative to adiabatic for 
an external wind (Chevalier & Li 2000). 

occurs in the Thomson regime for T > 10 — 3 day, i.e. the scattering cross-section 



With the aid of equations (|T^), (|20|), and (Q), it can be shown 
that the electrons are radiative if 



, 3(5 - s) (m e /m p ) 2 4.2 x 10 16 5 
nrT > 



32(F r + l) a e e e e B 



£&£b Y r + 1 



(37) 



which, with the further use of equations ( |l4| ) - (|l8|), leads to 
the conclusion that the electrons are radiative until the observer 
time T r given by 

(s = 0) T r = 0.025 E 53 n tfi (Y r + 1) 2 £ 2 _ x e\ _ 2 day , 

(38) 

(s = 2) T r = 0.23 A*(Y r + l)s e -xe B -2 day . (39) 

3.2.3. Adiabatic Electrons 

For 7, <C 7c equations ( j29| ) and ( |32| ) give the Compton pa- 
rameter 



Y n 



T e X 



7 r 1 7 c 3 - p 2<p<3 
7 2 3 < p 



(40) 



Case 1: 2 < p < 3. By substituting equations (|l9|) and ( J26| ) 
in equation (|40[), one obtains 

r a (F a + l) 3 ^ - T P (T) = c s {p)el- l ^-\nT 2 r w Y- 2 , 

(41) 

where log c s (p) — (3— p) log(5 — s) — log(3— s) + lAp— 3.7. 
The Compton parameter can be obtained by solving numeri- 
cally the above equation. For analytical purposes, one can ap- 
proximate Y a = T p for T p < 1, in which case the IC losses are 
less important, and 



Y„ 



■r-4-p 

■r p 



for T v > 1 



(42) 



In the latter case the IC scatterings affect the electron cooling. 

Note that the quantity nT 2 r in equation ( |4~ij ) decreases with 
time. Thus, for eb < £~b, the Compton parameter Y a is above 
unity until a time TL which can be determined by substituting 
equations @ - @ in @: 



(s = 0) T y = 10— £ 53 ™*,o£ 



1 r, 3 — » 



o £ ,-. i fc s,-3 2 day 



( S = 2) T y = 10— 



4 p p " 2 p 

/1 * e e,-l £ B 



3-P 
p-2 



day 



(43) 



(44) 



Thus, for eb < £~b and T r < T < T y , the Compton param- 
eter determines the evolution of the cooling break frequency 
(eqs. [|| and [||): 



For eb < £~b the Compton parameter is above unity until 

(s = 0) T y = 0.11 £ 5 3n*, <_i day , (48) 



(s = 2) Tj, = 0.87 A*£* _ x day 



(49) 



For eb < £~b and T, r < T < T y , the evolution of the cooling 
break frequency is 

(S = 0) V c = 1.1 X 10 17 ^53 3/2 ^K-l£B 3 -3 r i /2 Hz . 

(50) 

(,s = 2) v c = 1.5 x 10 16 eU'A^e-^Eb 3 ^ 2 Hz . 

(51) 

Note that in this regime v c increases with time in both models. 



3.3. Absorption Break 

The synchrotron self-absorption frequency v a can be calcu- 
lated with the aid of equation (6.50) from Rybicki & Light- 
man (1979). With the notations 7 P = min(7j,7 c ), v p = 
min(yi, v c ), and vq — max(^, v c ) it can be shown that optical 
thickness to synchrotron self-absorption can be approximated 
by 



(vjvp) 



-5/3 



V < u„ 



Bl b p I (vIvpY^' 2 V p<v<v a 



, (52) 



where £ = (3 — s) _1 nr is the remnant electron column den- 
sity, g = 2 for radiative electrons (j c < 7,), q = p for adiabatic 
electrons (74 < 7 C ). 

3.3.1. Radiative Electrons 



Equations (]20[), ([26]), and (52) give the optical thickness at 
the cooling break frequency 



= ^5_ enr = 1.1 xlQ-» (y^^/y^e 



3- Sj B 7 f (3-s)(5-s) 5 
For s = equations @, (|B]), and (||) lead to 

(S = 0) T C = 0.11 ^ 5 2 3^o + 1) 5 4 /2 



18 • 

(53) 



(54) 



Vr ( = 0) 10 15+ 



2.5p-5.5 

4 — p 



^53 ' t *,0 fc e,-l fc 



(s = 2) io 15+ " 



o3 



A~ 4 E 



2(P-1) = ." 
e,-l 



t B,-3- f <i 



Hz , 

(45) 

p 

Hz . 
(46) 



Note that for a homogeneous medium (s = 0) and 8/3 < p < 
3, the cooling break frequency increases with time, unlike the 
decreasing behavior it has for T < T r . 

Case 2: p > 3. This case is treated here for completeness, 
as there are no afterglows for which such a steep electron index 
has been found. Equations (|l9|), (f30|), and ( po| ) lead to 



-e e nr n 8 



(47) 



For t c < 1 the optical thickness to synchrotron self-absorption 

u 3/5 

is unity at v a given by v a = v c t c ' : 



s = 0) ^ a = 6.5 x 10 9 4 3 /10 n" /10 (^ + 1)4' -ilT" Hz . 



6/5 



-1/2 



For s = 2 equations (|TJ) - ([|]), and @ give 



(55) 



( S = 2) r c = 0.44 E 5 - 3/2 4l(r r + l) 5 £ ^ /2 T- 7/2 . (56) 



For T > T a we have r c < 1 and 2/ a < 

(s = 2) ^ a = 1.4x 10 12 -B 53 2/5 Ai 1/5 (y r +l) e | /5 _ 2 T ( 7 8/ 2 5 Hz 

(57) 
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3.3.2. Adiabatic Electrons 



The optical thickness to synchrotron self-absorption at the 
injection break can be found using equations (|l9|), (|20|), and 

& 



2.9 x 10" 



,1/2 



ri8 



1 3-sBt? 3-s ej 
If the external medium is homogeneous 

= 0) n = 1.3x10 
( S = 0) 



1/2 
-I s B 



r e 



p — 1/ 2 „ _ — O — >-l"rr\ 



5 -1/2^5/2 
-l fc B,-2-'d 



^ a = 2.6 x 10 9 S 5 1 3 /5 n^ 5 e -^ 1 e^ /5 _ 2 Hz 



which is time-independent. 
For the wind model 



6.8 x 10~ 7 E, 



-3/2 
53 



42.-5 -l/2 r 3/2 



(58) 

(59) 
(60) 

(61) 



3.7 x 10 9 E^ATe-Ue^T-^ Hz 



( S -2) r, 
(s = 2) ^ a 

(62) 

Note that, in general, the absorption frequency decreases 
faster for a remnant interacting with a wind than for one run- 
ning into a homogeneous external medium. 

4. ANALYTICAL LIGHT-CURVES 

If the effects arising from the remnant spherical shape (see 
Appendix) are ignored, than the observed flux peaks at v p = 
min(i/i, v c ), where it has a value 



F u 



tTBN c 



(63) 



Here 



4irD 2 m e c^ 

bp is a factor calculated by Wijers & Galama (1999), 



which we shall set p = 0.63, D = (1 + z)-^ 2 Di{z) with 
Di the luminosity distance, and N e = m{r) /nip is the number 
of electrons in the remnant. Equations (|l9|), (pCh, and (|53|) give 



F u 



57 



n-2 l/2p2 3/2 3 
iJ 2S, t B 1 n '18 



mjy 



(64) 



The afterglow emission at any given frequency and time 
can be calculated using the synchrotron spectrum for the elec- 
tron distributions given in equations (|T]) and ( f32"| ) (e.g. Sari et 
al. 1998): 



F v = F Vr 



{V/Va) 2 {v a /Vc) l t* 

(ylv c )- 112 



V < V a (1) 

V a < V < V c (2) 

v c < v < Vi (3) 

Vi < v (4) 



(65) 



for T < T r , assuming that v a < v c , and 



F v = F^ 



(iz/i/i)-^- 1 )/ 2 



v <v a (5) 
v a < v < vi (6) 
Vi < v < v c (7) 
v c <v (8) 
(66) 

for T > T r , assuming that v a < v^. 

In Figures 1 and 2 the plane T — es is divided into several 
regions which are labeled as in equations ( |65| ) - (j66|), according 
to the ordering of the observing frequency v and of the three 



break frequencies v a , Vi, and v c . The observed fluxes in each 
case are given in the Appendix, and a set of multi- wavelength 
light-curves is shown in Figure 3. The results shown in Figures 
1^1 have been obtained using equations which are valid in any 
relativistic regime, such as equations (Q), (19), (p0|), and (63). 



4. 1 . Inverse Compton Emission 

The IC emission can be easily calculated by using the above 
equations for the synchrotron spectrum and Compton parame- 
ter. The up-scattered spectrum peaks at z/, c ~ 7^ v c if electrons 
are radiative and at vi c ~ -ffvi if electrons are adiabatic. It can 
be shown that, for any electron radiative regime, the flux of the 
up-scattered emission at this frequency is 



(67) 



where r e and F„ p are given by equations ( |3C| ) and (|64|), respec- 
tively. For the up-scattering of synchrotron photons above v a 
and assuming that Y a < 1, the resulting IC light-curves have 
the following behaviors: 



s = 

T < T r 



T > T r 



s = 2 

T < T r 



s = 2 

T > T r 



oc 



p( ic ) 



? (ic) 



jrl/3 




v £ 7>c 




T l/8 




l 2 c v c <v< 




y-(9p- 


10)/8 


l 2 v l < V 


(68) 


T 1 




v < jfvi 




j>— (9p- 


U)/8 


itvi <v< 


7>c 


j>-(9p- 


10)/8 


llv c £ V 


(69) 


p-5/3 




$ 7>c 




rpO 


7 C 2 


v c <v< 7? 


^ , 


rp~(p- 


1] if 




(70) 


T -l/3 


V 






rp-p 


7 2 


Vi<v<-i 2 c 




JT-O- 


1} 7 C 2 


v c < V 





(71) 

For external media that are not denser than assumed so far, 
the IC emission is weaker than synchrotron, even in soft X-rays. 
As shown in the upper left panel of Figure 3, for n* <; 10 cm~ 3 
and A* <; 1 the up-scattered radiation can dominate the syn- 
chrotron emission at times T > 10 _1 day, diminishing the de- 
cay rate of the X-ray emission. This is due to that the flux at 
the IC peak (eq. J67|]) depends strongly on the external medium 

density: Fu^ oc nj^ for s = and Ful^ oc ^4*^ 2 for s = 2. 
The afterglow flattening is strongly dependent on the observing 
frequency, being absent in the optical and below. 

5. CONCLUSIONS 



Using the analytical results given in equations ( |B1[ ) - ( JC9| ), 
the afterglow light-curve can be calculated at any frequency and 
at observing times up to the onset of the non-relativistic phase. 
As illustrated in Figure 3, the largest differences between the af- 
terglow emission in the two models for the external medium is 
seen at low frequencies (lower panels). However, the scintilla- 
tion due to the local interstellar medium (Goodman 1997), may 
hamper the use of the radio light-curves to identify the type of 
external medium and geometry of the ejecta (Frail et al. 2000). 

Figure 4 shows that, for various model parameters, the rate 
of change of the afterglow emission at v ~ 10 12 Hz and at 
early times (when the jet effects are negligible, provided that 
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the jet is initially wider than a few degrees) exhibits a strong 
dependence on the type of external medium. If the external 
medium is homogeneous the sub-millimeter afterglow should 
rise slowly at times between ~ 1 hour and ~ 1 day, while for a 
pre-ejected wind the emission should fall off steeply, followed 
by a plateau^j Therefore observations made at sub-millimeter 
frequencies with the SCUBA (James Clerk Maxwell Telescope) 
or with the MAMBO (IRAM Telescope) instruments would be 
very powerful in determining if the medium which the remnant 
runs into is homogeneous or follows a r~ 2 law. 

We note that if turbulence in the shocked fluid does not lead 
to a significant mixing, then the inhomogeneous electron dis- 
tribution will alter the afterglow spectrum below the absorption 
frequency v a as described by Granot, Piran & Sari (2000). The 
result is that the afterglow emission at v < v a rises more slowly 
than calculated here. For instance, the T 2 rise exhibited by the 
v = 10 12 Hz light-curves shown in Figure 4 for the wind model 
at early times becomes a T 1 rise. Nevertheless, the basic differ- 
ence mentioned above between the temporal behaviors of the 
sub-millimeter light-curves at 10 -2 day <T<1 day remains 
unchanged. 

The IC losses alter the evolution of the cooling break v c if 
the electrons injected with minimal energy are adiabatic and if 
the Compton parameter is above unity (i.e. the magnetic field 
parameter is weaker than that given in eq. [36|). In the case of 
a homogeneous external medium, the cooling break frequency 
decreases as T -1 / 2 if the electrons are radiative. When the 

3p-8 

electrons become adiabatic, this break evolves as T 8 - 2 ? for 
p < 3 and increases as T 1 / 2 for p > 3. For an external wind 



the change is from v c oc T 1 / 2 to v c oc for p < 3, and 

to v c oc T 5 / 2 for p > 3. Consequently the power-law decay of 
the afterglow emission at frequencies above the cooling break 
flattens by up to 1/2 if the external medium is homogeneous 
and by up to 1 if the medium is a wind. For an electron index 
p < 3, the flattening is mild and likely to be seen only in the 
optical emission from a remnant interacting with a pre-ejected 
wind. 

The IC emission itself is generally weaker than the syn- 
chrotron emission. Nevertheless, if the external medium is suf- 
ficiently dense (i.e. n* ^ 10 cm~ 3 or A* ?Z few), a flattening 
of the soft X-ray light-curve should be seen few hours after the 
main event, at fluxes well above the threshold of BeppoSAX 
(see Figure 3, upper left panel). The flattening of the after- 
glow emission due to the up-scattered radiation is a chromatic 
feature, appearing only at high frequencies, and its strength is 
moderately dependent on the remaining model parameters. 

Finally, another possible signature of the interaction with the 
wind ejected by a Wolf-Rayet star should be found during the 
GRB emission in the form of smooth pulses of increasing du- 
ration. Such pulses are generated in internal shocks when the 
decelerating outermost shell is hit from behind by shells ejected 
at later times. This phenomenon is more likely to be seen in the 
wind model, for which the deceleration radius is smaller than 
for a homogeneous medium. 

AP acknowledges the support received from Princeton Uni- 
versity through the Lyman Spitzer, Jr. Fellowship. We thank 
Bohdan Paczyriski for useful discussions. 



APPENDIX 

The observer time calculated with equation ( p"2| ) is the time when photons emitted by the shocked gas moving precisely toward 
the observer arrive at detector. Photons emitted by the fluid moving at a non-zero angle relative to the direction toward the observer 
are less boosted by the relativistic motion of the source and arrive at the observer later. The effect of the remnant spherical curvature 
on the afterglow emission can be reduced to a correction factor that must be applied to the light-curve obtained analytically using 
the remnant parameters at laboratory frame time t related to the observer time T through equation (|l2|). This correction factor is 
dependent on the observing frequency, as described below. 



CORRECTION FACTORS FOR ANALYTICAL LIGHT-CURVES 

The flux received by the observer at time T can be calculated by integrating the remnant emission over the equal arrival-time surface 
(Panaitescu & Meszaros 1998). This surface is defined by cT = ct — r cos 6, where 9 is the azimuthal angle measured relative to 
the observer's line of sight toward the remnant center. Using the equation for dynamics of relativistic remnants (r oc r~( 3 ~ s )/ 2 ), the 
integral can be written as 

A{cT) -2 f R(T) (p^ )er 4- S 



F ^ T) -w^L Hrnr*' (A1) 

where x = (3 — s)(r / R) 4 ~ s , and R is the radius for which photons emitted at 9 = arrive at T, i.e. the maximal ra- 



dius on the equal arrival-time surface. In equation (Al) (P^/) e is the co-moving frame power per electron, taking into ac- 
count the spectrum of the emission. For the synchrotron spectra given in equations ( j(55[ ) and (|6g), one can write generically 

(P' u i)e = ^/i<t>p{e i B / m e c 2 )(v' I v' n ) aa (is' /v' i ) ai (f' / y' c ) ae , where v' is the co-moving frame observing frequency. The integral 



in equation ( Al ) is determined by the the values at of the integrand at r % R. 



If the effect of the remnant geometrical curvature on the photon arrival-time were ignored, then the observed flux would be 

Fl°\R) = [(Pl,)eTN e ] R{T) , (A2) 



which can be calculated based on equations (p3|), ( p5[ ) and (|56|). The r-dependent quantities in equation ( |A1[ ) can be expressed in 

2 The light-curves presented in the upper left panel of Figure 4 show that this criterion for determining the type of external medium fails only if the particle density 
of the homogeneous medium exceeds ~ 10 cm -3 . In this case the X-ray emission may help to distinguish between the two models of external medium, as the 
absence of a flattening of the high energy emission is compatible only with a pre-ejected wind with A, less than a few. 



terms of their values at r = R, so that the flux F v can be written as times a correction factor 



K = 2(4- s) 2 - iQ « +a > +a ^ f u fis) [1 + (3 - s)u 4 ~ s ] 

Jo 



4— si a«+Oi+« c -2 



(A3) 



where f(s) = 7 — 2.5s — 2a a (q — 3s) + c^(2 — 0.5s) — a c (2 + 0.5s), with q = 1, 2 for radiative, radiative electrons, respectively. 
The X factor depends on i/ and also onp if v > z/j. Table 1 gives its values for various cases. 



HOMOGENEOUS EXTERNAL MEDIUM (S = 0) 



By substituting the equations for the break frequencies and equation ( p4| ) in equations ( |65| ) and (66), and taking into account the 
above correction factors for the remnant curvature, the following fluxes are obtained 



(i) 



F v v =' 0.3 Dzi(Y r + l)- X n-lez]_ x v^T^ mJy , 
F v ® 10 D;i(Y r + l)^EU 6 n^e B ,_ 2 ^»T^ a mJy , 
F v ® 40 ^(Y, + l)-^5 3 3 /4 ^^ ^f^ 4 mJy , 
F v i i a-i-o.«P D -2 (Fr + el:\ef_ x v^T^ mJy 

F„ ® 30 D^E^n-y^,-! vl. 7 T^ mJy , 
F {6) I D- 2 E 5/6 n 1/2 £- 2/3 £ 1/3 ^/ 3 T 1/2 mJv 

r " ~ 1 -^28 ^53 ' t *,O fc e : -l t B,-i V IA.6 1 d,-2 IUJ y ' 

p+i _ p — 1 



^ (X) 10 2.i-i. 3p n-aR-T^V^-i^ ^t:^- 1 

(8) 



p + 2 p-2 _p 

) — 2 771 4 c-P - 1 c 4 „ 2 T 

'28 -^53 fc e,-l fc _B,-2 "l4.6 J d 



mJy 



F v ^ IT^^ 



1(12-/) A, 
-^53 'Vo 



-2) (p-l)(3-p) i(-p 2 +2p+4) 
fc e,-l t B,-4 



_ p _ il£ 4- l 
"l7.5 J d,-l 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 

mJy (2 < p < 3) . (B9) 



The case given in equation (B9) and labeled (8a) corresponds to the same frequency ordering as for case (8), but the cooling break 
v c evolution (eq. [p3[]) is determined by the IC losses, i.e. T r < T < T y and Y a > 1. 



WIND EXTERNAL MEDIUM (S = 2) 
Following the same exercise as above and using the relevant equations, the following results can be obtained for the wind model: 



F v ( = 3 0.03 D^ s 2 (Y r + ly'EssA-^- 1 ^ e^-2*M-i mJ y . 
F v { B 70 D- 2 2 {Y r + l)^EU 3 A 5 J S e B ,^ v\£t£!* mJy , 
F v ( =' 0.07 D^EsaA-he,-! ^l-i mJ y . 

- 9 ^28 ^53 A * £ e -1 £ S,-3 ^12 J d mJ y J 



F„ ( = } lO 2 - 3 " 1 -^ D^Eli' A^-\e^. v^T^ ' mJy 



Sa) . 1.9p^-8.6p + 5.4 



= 10 *=5- 



^28 -^53 



— f 4 -(p-2) _(p-l)(3-p) = .i(-p 2 +2p+4) 

-^i* c„ _i 



( P) mJy (2 < p < 3) 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 

(C6) 



For the remainder of the cases it can be shown that the fluxes for the wind model differ from those for obtained in the s = model 
only by a constant factor: 



F„(s = 2) (3) 06Q F r (s = 0) 



F v (a = 0) 



Y r (s = 2) 



F u (s = 2) (4) /l7V /4 y r (s = 0) 
~ \72j Y r (s = 2) 



F v {s = 0) 

F v {s = 2) (8) 
F v (s = 0) 



1.35 



-) 

72 J 



(C7) 
(C8) 
(C9) 
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This means that observations made at frequencies above the cooling break are unable to distinguish between the two models for the 
external medium, with the exception of the case where the electron cooling is dominated by IC losses (case 8a). 

The reference frequency in the equations above was chosen v = A x 10 14 Hz whenever this optical frequency falls in one of the 
cases labeled in equations (|65|)-(|66|), for the model parameters given in Figures 1 and 2. The observer time was scaled to a value at 
which each case is more likely to occur. 

The above light-curves must be corrected to account for cosmological effects. This is achieved by replacing T with T/(l + z) and 
v with (1 + z)v, where z is the burst redshift. 



TABLE 1. 



Values of the factor K given in equation (|A3|), representing the correction that should be applied to the 

ANALYTICAL LIGHT-CURVES TO ACCOUNT FOR THE EFFECTS DUE TO THE REMNANT GEOMETRICAL CURVATURE. 



s 


(1) (2) 


(3) 


(4)* 


(5) 


(6) 


(7)* 


(8)* 





0.50 0.55 


0.59 


1.37 


1.00 


0.42 


1.34 


1.02 


2 


0.11 1.76 


0.74 


1.19 


0.12 


0.90 


2.25 


0.82 


case 


K depends 


on p. 


The 


values 


given 


here 


are for 
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FIG . 1 . — The afterglow brightness at a given frequency (indicated in each panel) is dependent on the relationship between this frequency and those of the injection 
(ui), cooling (u c ), and absorption (u a ) breaks. Different combinations iiLthese_four frequencies are indicated in the time — eb plane for the homogeneous external 
medium model (s = 0), according to the labeling given in equations M) - The label 8a indicates the case when v-i < u c < u and Y a > 1, i.e. the cooling 
break v c evolution is determined by the IC losses. Other model parameters are E = 10 53 erg (in 47r sr), n* = 1 cm -3 , e e = 0.1, and p = 2.5. The upper two 
panels also indicate the time when the Compton parameter Y a for adiabatic electrons becomes unity. 
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FlG. 3. — Numerical light-curves for D = 10 28 cm, E = 10 53 erg (in 4n sr), e e = 0.1, £b = 10~ 3 , p = 2.5, for a homogeneous external medium (dashed 
lines) with n» = 1 cm'jL^anijLa.pre-ejected wind (dotted lines) with A* = 1, at various observing frequencies (given in each panel). The numbers indicate the cases 
identified in equations (p5j) - Q6q). The afterglow emission is calculated using equations that are valid in any relativistic regime. The effect of the remnant geometrical 
curvature on the photon arrival time and relativistic boosting are taken into account. 

Note that the largest differences between the two models are shown by the low frequency light-curves at observer times below ~ 1 day. 

Also shown in the upper left panel are two soft X-ray light-curves for denser external media: n» = 10 2 cm -3 for s = (thin dashed curve) and A, = 10 for 
s = 2 (thin dotted line), in which cases the afterglow emission at T 10 day is dominated by IC scatterings, and exhibits a substantial flattening. Such a feature 
is chromatic and cannot be seen at lower observing frequencies. The flattening of the X-ray emission in the case of denser media is present for other plausible sets of 
model parameters, as long as they lead to fluxes (at the time when the flattening occurs) that are detectable with current instruments. 
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FIG. 4. — Comparison between the afterglow light-curves at observing frequency v = 10 Hz for homogeneous external media (dashed lines) and pre-ejected 
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